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' We use the non-perturbative renormalization group to clarify some features of 

perturbation theory in thermal field theory. For the specific case of the scalar field 

O; 

theory with O(N) symmetry, we solve the flow equations within the local poten- 
tial approximation. This approximation reproduces the perturbative results for the 
screening mass and the pressure up to order g 3 , and starts to differ at order g 4 . The 
JL method allows a smooth extrapolation to the regime where the coupling is not small, 

very similar to that obtained from a simple self-consistent approximation. 
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PACS numbers: 05.10.Cc, ll.10.Wx, 11.15.Tk 

I. INTRODUCTION 

Perturbative calculations at high temperature suffer from bad convergent behavior. Even 
the thermodynamical quantities, such as the pressure or the entropy density, which are 
dominated by short wavelength degrees of freedom - for which one could expect perturbation 
theory to be reasonably good - are not well described by a strict expansion in powers of the 
coupling constant. This problem has been particularly studied in quantum chromodynamics, 
in the context of the quark-gluon plasma (for a recent review, see l|). In this case, the 
perturbative expansion of the thermodynamic potential is known up to order g 6 log g, where 
g is the gauge coupling (the g 6 term requires non-perturbative methods {(]]). 
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However, the problem is not specific to QCD at high temperature: Similar poor convergence 
behavior appears also in the simpler scalar field theory [7], and has also been observed in 
the case of large- N 4 theory jg| . 

Reorganizations of the perturbative expansion, based on various arguments, have been 
proposed in order to extend the usefulness of weak coupling calculations to regimes of not too 
small coupling. These involve "screened perturbation theory" |9|], "hard thermal loop pertur- 

, k], [ijl , and the approach put forward in Refs. 



bation theory" [10, [1J 
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171 ]. based 



on an expansion of the thermodynamical potential 
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in terms of dressed propagators which 



Q, 3- The approach of 1 
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involve two-particle irreducible (2PI) diagrams 
exploits a nonperturbative expression for the entropy density that can be obtained from a 
^-derivable two-loop approximation 211 ] . Here the emphasis is on a physical picture involv- 
ing quasiparticles whose residual interactions are assumed to be weak after the bulk of the 
interaction effects have been incorporated in the spectral properties of these quasiparticles. 
It has been shown in particular that the entropy density of the quark gluon plasma can be 
well understood within such a scheme down to temperatures T > 3T C where the coupling 
can be as large as g ~ 2 |17j ]. As an alternative to resummations, effective field theories 
have also been used [4] . Among those, dimensional reduction, which emphasizes the role of 
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231 ] . An interesting 



the zero Matsubara frequency, stands out as a very powerful one 
feature of the perturbative calculations when organized through the dimensionally reduced 
effective theory is that the large scale dependence of strict perturbation theory is consid- 
erably reduced when the effective parameters are not subsequently expanded out 
When this is done, the predictions of dimensional reduction become similar to those of the 
2PI resummation. 

In fact, the origin of the difficulties of thermal perturbation theory is now well understood. 
What complicates the situation is the fact that the coupling constant alone does not control 
the magnitude of corrections to the free theory: Thermal fluctuations of various wavelengths 
play also an essential role. In the weak coupling regime, the thermal system is characterized 
by a hierarchy of scales: Most particles have momenta of the order of the temperature. At 
high temperature, these "hard" degrees of freedom dominate the thermodynamics and their 
interactions are accurately described by ordinary perturbation theory. However, perturba- 
tive corrections to the thermodynamical functions involve also "soft" degrees of freedom, 
whose momenta are typically of order gT. These soft modes are non-perturbatively related 
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by their coupling to the hard modes. Such corrections are easily handled by resummations 
or effective theories. Furthermore, the soft modes also interact among themselves. Although 
perturbative, these interactions generate corrections which take the form of an expansion 
in powers of g rather than of g 2 as is the case for ordinary perturbation theory: as a re- 
sult, perturbation theory becomes less accurate. In QCD, there exists a further scale, g 2 T 
for which perturbation theory stops to make sense: at this scale, the self interactions of 
the modes become comparable to their kinetic energies, invalidating any expansion around 
free particle motion. The problems associated with this "ultra-soft" scale are specific to 
non-abelian gauge theories; they will not be discussed in this paper. 

There are indications, already alluded to above, that the structure identified in weak 
coupling, which is based on the hierarchy of scales that we have just discussed, seem to sur- 
vive, after appropriate resummation, or appropriate use of effective theories, even in regimes 
of strong coupling where the arguments used to identify the structure become unjustified. 
That is, the resummations that have been motivated by analyzing the weak coupling regime, 
provide a smooth extrapolation into the regime of strong coupling where a strict expansion 



in powers of the coupling does not make sense. The purpose of this paper is to shec 
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this issue by using the non-perturbative renormalization group (NPRG) 

There is some analogy between the effective field theory approach and the NPRG: in 
effective field theory one integrates out degrees of freedom above some cut-off; in the renor- 
malization group this integration is done smoothly. In a sense, the renormalization group 
builds up a continuous tower of effective theories that lie infinitely close to each other and 
are labeled by a momentum cut-off scale k. These effective theories are related by a renor- 
malization group flow equation. The picture remains essentially the same for any value of 
the coupling. 

The NPRG has been applied, in various incarnations (it is also called exact, or functional) 



to a variety of problems in condensed matter 



and nuclear physics 



temperature [27 
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(for reviews see [39 
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r[30 
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311 ]. in particle [32 
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). It has been applied to problems at finite 
461 ] . and the general behaviors that we shall report in this 



paper have been known already for some time. However, the main focus of previous studies 
has been the description of the phase transitions, rather than the specific problem that we 
want to address here. 

The outline of this paper is as follows. In Sec. [TT], we present a brief review of available 
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results concerning perturbation theory for a scalar field with O(N) symmetry. We also 
discuss a simple self-consistent ^-derivable approximation that becomes exact in the large 
N limit. Section [Till gives a brief introduction to the non-perturbative renormalization group 
and the local potential approximation. Specific features of finite temperature calculations 
are recalled. In particular we introduce a regulator that we found particularly convenient 
for such calculations. In Sec. [IV] we integrate the flow equations numerically and discuss 
the results obtained. Finally, in Section [V] we perform a perturbative analysis of the flow 
equations and comment the results. Note that a preliminary account of this work was 
presented in Refs. |47|, 



II. PERTURBATION THEORY AND THE NEED FOR RESUMMATION 

In this section we briefly review existing results of perturbation theory for the thermody- 
namics of the scalar field. We also recall how a simple self-consistent approximation based 
on the lowest oder 2PI diagram, can be used to include non perturbative effects that allow 
for a smooth extrapolation to strong coupling. This will be useful later in our discussion, as 
it turns out that the results of this self-consistent approximation are close to those of the 
non-perturbative RG within the local potential approximation. 



A. Brief review of results from perturbation theory 

We consider a scalar 4 theory in d dimensions, defined by the O(N) symmetric La- 
grangian 

C(x) = \ [d V (x)f - l -m\\x)-g> [p\x)]\ (2.1) 

where the field (p(x) has N real components <fi{x), with i = 1, ...,N. The temperature T 
{(3 = 1/T) enters through the action 

S = dr j d d ~ l xC(x), (2.2) 

and the periodicity of the fields in imaginary time, (p(/3) = <f(0)- We shall fix the bare 
parameters of the action at T = 0, in order to have a vanishing renormalized mass and a 
given value of the renormalized coupling constant (the last point is discussed in Sec. Ill Bl) . 



Figure 1: Diagrams that contribute to the screening mass m 2 D at order g 2 (left) and at order g 3 
(right). 



Then we turn on the temperature at a given fixed value and we study the finite temperature 
predictions of the model for various values of the zero-temperature coupling. 

There are two physical quantities that we shall discuss in this paper: the screening mass 1 , 
and the pressure. Up to order g A , the screening mass reads 

ft 



23|]: 



In 



2ttT 



3 7T 



2.415 
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(2.3) 



where p, 2 = 4ne~ J fi 2 with /x 2 the scale introduced in dimensional regularization, 7 the Euler 
constant, and here g 2 = g 2 (fi) is the MS coupling constant at the renormalization scale ft. 
Similarly, the weak-coupling expansion of the pressure has been computed (in the massless 
case) to order g 5 , and reads [3, 
T 15 

P = Po 



23. 



49) 



--ff) 2 +Y(r 
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(2.4) 



where P = (tt 2 /90)T 4 is the pressure of an ideal gas of free massless bosons. 

The lack of convergence of the weak-coupling expansion of both the screening mass and 
the pressure, manifests itself by the fact that, unless the coupling is very small, g ^ .1, 
the successive corrections do not decrease in magnitude, and the dependence of the renor- 
malization scale gets larger and larger (see e.g. This is in contrast to what one would 



1 Note that the screening mass, defined as the pole of the static propagator (for complex wave numbers), 
differs from the quasiparticle mass, defined as the pole of propagator at vanishing three-momentum. Up 
to order g 3 , the two quantities coincide in scalar theory but start to differ at order g 4 . Furthermore, 
within the approximations that we shall consider in this paper, where the momentum dependence of the 
self-energy is neglected, these two masses are the same, and coincide also with the second derivative of 
the effective potential with respect to the field. 
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expect in a well behaved perturbative expansion where the explicit scale dependence cancels 
against that of the running coupling, to within terms of higher order than those retained in 
the calculation. Such cancellations do indeed occur in Eqs. ( 12.31) and (12 Ah . To see that, it 
is enough to use the one-loop j3 function: 

"f^ 4 + <^ < 2 - 5 > 

However, since the successive orders do not get smaller and smaller, the cancellation of the 
scale dependence remains only formal. 

It is instructive to recall the origin of the first two terms of the expansions of the screening 
mass and the pressure (which can be obtained from the diagrams depicted in Figs. [1] and 
[2]), because they illustrate well some aspects of the physics that we are discussing. We shall 
work out only the calculation of the mass; that of the pressure is similar. 

The first two terms in Eq. ( 12.31) are obtained from the Feynman diagrams displayed in 
Fig. [T] which we calculate as (using an ultraviolet (UV) cut-off rather than dimensional 
regularization for the vacuum part) 



= lW A ^±|^ + ^, (2,) 

with u 2 = q 2 +m 2 . In the second line of ( I2.6P the sum over Matsubara frequencies u n = 2irnT 
has been converted into an integral over the distribution function n(u) = [e^l T — 1) _1 (see 
Appendix A), and a UV regulator A introduced. The counterterm 5m 2 , of order g 2 , cancels 
the ultraviolet divergence (as A — > oo) of the vacuum integral and its finite part is chosen 
so that the renormalized mass vanishes in the vacuum. 

The order g 2 contribution is a genuine perturbative correction, dominated by the hard 
degrees of freedom. It can be estimated by neglecting the mass in u q : 

,2 



m 



D 



% ^ d qq n{q) = g 2 T 2 . (2.7) 
71 Jo 



The g 3 term however is not, strictly speaking, a perturbative correction: the odd power is 
the result of an infinite resummation. What is at work here is precisely the coupling between 
soft modes and hard ones that we discussed in the introduction: when the momentum 
running in the loop of the left diagram of Fig. [I] is soft, i.e., of order gT, the correction to 



Figure 2: Diagrams contributing to the pressure at order g 2 (left diagram) and g 3 (right diagram). 



the propagator due to the hard fluctuations cannot be ignored since these hard fluctuations 
contribute a mass also of order gT. Thus one needs to keep the thermal mass in the 
propagator. Starting from (12. 6p and subtracting the 0(g 2 ) contribution that we have just 
calculated, we obtain 

lg 7T 2 Jo \ u q q J 

= --g 3 T 2 . (2.8) 

7T 

In the second line we have used the approximated form of the statistical factors n(u) ~ T/u, 
appropriate since the integral is dominated by soft momenta q <C T. (One would have 
obtained the same second line of (12.81) by starting from the expression in the first line of 
(12. 6p and keeping only the u n = contribution.) 

The type of resummation involved to get the g 3 term can be turned into a fully self- 
consistent approximation which extrapolates smoothly to strong coupling. But before we 
do that, let us add a few remarks about the running coupling, and in particular about what 
is meant by strong coupling in the present discussion. 



B. Remarks about the coupling 

The running of the coupling, to order one-loop, is governed by the /5-function (12.51) . 
Assuming an ultraviolet cut-off A, where the value of the coupling is g& and integrating 
Eq. (12.51) . one gets: 

1 1 9 , A , , 

il^ + ^ ,0g ? (2 ' 9) 
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Figure 3: Diagrams contributing to the four-point function at order g 2 (left diagram) and g A (right 
diagram). 

which shows that — > when [l — > 0. The coupling at scale /i = is therefore not a suitable 
quantity to characterize the theory. In the following, we therefore adopt the usual practice 
of fixing the coupling at the scale 2ttT, which we shall simply denote by g = g(p, — 2ttT) 
(unless ambiguities may arise). Equation (12.91) can then be rewritten as 

1 1 9 2nT 

Now, for any finite value of g there is a scale fi = Al (the "Landau pole") at which the 
coupling diverges: 



A L = 27rTexp — . (2.11) 



2tt_ 2 

If g <C 1 then Al ^> A and the physics is not affected by Al. However, as g grows, 
Al decreases and eventually becomes of same order as A. In order to avoid unphysical 
behaviors, we require A L > A, which puts a constraint on the largest admissible values of g: 

2n 2 

It is easy to verify that this bound corresponds to an infinite value of g\. To increase the 
maximum value of g, one could increase the temperature; however if 2irT becomes too close 
to A, the results become sensitive to the value of the cut-off (see the section on numerical 
results for further discussion of this point). Note that this estimate of the upper bound for 
g is modified by non-perturbative effects, as we shall see in Sec. IIVI 

The coupling constant may be viewed as the value of the scattering amplitude A for van- 
ishing external momenta. While, as we have just seen, after renormalization this scattering 
amplitude vanishes, this is not so at finite temperature. Up to order g A , A can be calculated 
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from the diagrams in Fig. [31 We obtain 

A = g 2 + 36g A TY [ -pL ( \ r V . (2.13) 

»J ( 27r ) Wn + q 2 + ™ 2 D J 

We can calculate the sum-integral using the formulae given in Appendix A or, more di- 
rectly, noticing that the thermal contribution is infrared dominated, by keeping only the 
contribution of u n = in (12.131) . We obtain 



., 1% 4 f°° 2 T 2 9g 3 
A^g l -^r\ qdq—=g 2 --^, 2.14 
7T 2 ./o 2ir 



where g 2 = g 2 (2irT). Thus the scattering amplitude contains a g 3 contribution, whose origin 
is the same as that of the g 3 contributions in the mass or the pressure. 

C. Self-consistent 2PI resummation 

As we have seen earlier, resummations are required in order to deal with the IR aspects 
of thermal perturbation theory. We shall briefly recall here the equations that are obtained 
in the simple ^-derivable approximation where the single 2PI diagram that is considered is 
the left diagram in Fig. It is easily verified that this approximation is exact in the large 
N limit, in the sense that it realizes the resummation of all the diagrams that survive this 
limit. 

With this method, one obtains the following renormalized self-consistent gap equation 
for the mass m 

d 3 k n(e k ) 3^m 2 ( ^ m 2 



where the renormalized coupling is related to the bare one gs by (e = 4 — d): 



1 ^ 3 



By using Eq. (12.161) . one can check that the solution m 2 of Eq. (I2.15P is independent of \x 
(to all orders in g). It is easily seen that the /5-function that is deduced from Eq. (I2.16P 
is only one third of the one- loop /9-function (see Eq. (12.51) ). This is a peculiar feature of 



^-derivable approximations 50|: the particular resummation involved in the solution of the 



gap equation corresponds, for the scattering amplitude, to the iteration of the basic vertex 
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in a single channel out of three. The correct one-loop /^-function is recovered when one keeps 

n 

in the $- derivable approximation the skeleton of order g 4 [50j. 
One can also easily compute the pressure 

d 3 k _a £ m 2 f d 3 k n{ek) m 4 



p = - T J log(1 " e "*"» + T J Wf If + 12^ ■ < 2 - 17 ' 

One can verify explicitly, but this is obvious from a simple analysis of the diagrammatic 
content of the approximation, that this expression for the pressure is perturbatively correct to 
order g 3 . Besides, the complete, self-consistent results, as obtained by numerical evaluation 
of Eqs. (12.151) and (I2.17p . extrapolate smoothly to large values of g. We present the numerical 
results of the self-consistent mass and pressure in Sec. [IV] (in Figs. [13] and ]T5|) when we 
compare them to the results obtained in this paper by RG techniques. 



III. NON-PERTURBATIVE RENORMALIZATION GROUP 



The basic strategy of the NPRG [40|, ]5l|, [52| consists in adding to the classical action 
(12. 2D a regulator, conveniently defined in momentum space as 



&S K [<p] 



T 



E 



d d ~ l q 



(2n 



d _j<Pi(.q)R lt (q)<Pi(-q), (3.1) 



where q — (q — iu n ,q), with the Matsubara frequencies ui n = 2-rmT and R K is a cutoff 
function, depending on the continuous parameter k, whose specific form will be given shortly. 
The role of AS K is to suppress the fluctuations with momenta q < k, while leaving unaffected 
those with q > k. Thus typically R K (q) — > k 2 when g C k and R K (q) — > when q ^> k. 
There is a large freedom in the choice of R K (q), abundantly discussed in the literature 



53, 



54 



55 



56. 



57|, 
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60]. 



In this paper, we use a regulator which depends only on the spatial components of the 
momenta, q, but not on the energy variable go- It is chosen of the form 



^ h (q) 



(3.2) 



A regulator that preserves Euclidean invariance, and is of the form of Eq. (13. 2\\ with q — > q, 



could of course be used 



56j. However, such a regulator leads to several difficulties: since 



the regulator cuts off frequencies sharply, the contribution of single Matsubara frequencies 
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K = A 



H = 2vrT 



T = 




bare quantities, 
scheme dependent 



T > 



k = k = physical quantities 



Figure 4: Basic concept of how to apply the flow equation at finite temperature. The coupling g 
is extracted from the vacuum flow at the scale k = [i with /j, = 2ttT. 

become visible as we vary k, resulting in an oscillatory behavior of the flows of various quan- 
tities. Furthermore, because the fluctuations are not entirely suppressed at the microscopic 
scale A by such a regulator, initial conditions for quantities like the pressure need to be 
fine-tuned with temperature dependent counterterms. Although this can be done, it leads 
to unnecessary complications. 

The NPRG constructs a family of effective actions T K [<f)] (with the expectation value 
of the field), in which the magnitude of long wavelength fluctuations are controlled by the 
regulator R K . The effective action interpolates between the classical action obtained 

for n = A (with A the microscopic scale at which the fluctuations are essentially suppressed) 
and the full effective action obtained when all fluctuations are taken into account, that is, 



when k — > 0. One can write for T K [(j)] an exact flow equation [271 . l28l . |29| 

(2) 

where is the second derivative of T K with respect to 0. 

Conceptually, we will follow the strategy developed in Ref. 61| and apply the flow equa- 
tions in the way illustrated in Fig. HJ Starting with given physical parameters at k — 0, we 
integrate the flow equations up from k = to A thereby removing quantum fluctuations step 
by step in order to arrive at bare quantities at a chosen scale A. If A is chosen big enough, 
only renormalizable parameters survive and the system can be described by a simple set of 
bare parameters. Starting from these bare parameters we can follow the flow down from 
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k = A to 0, but this time with the temperature T turned on. The physical quantities are 
then obtained at k — 0. 

Note however that the coupling constant is fixed at the scale k = 27rT on the T = flow, 
as discussed above. As for the mass at scale A, it is adjusted so that the mass at k = 
vanishes. Note that this procedure induces a specific scheme dependence attached to the 
choice of the regulator. We should keep in mind this scheme dependence when comparing 
with results of perturbation theory or those of the 2PI resummation, that are obtained in 
another scheme. 

Computationally, it might be numerically hard to integrate the flow equations up from 
k = to A. But this is just equivalent to integrating them down, while carefully adjusting 
the bare parameters such that the flow will arrive at physically desired quantities at the end 
of the flow at k — 0. 



A. Local potential approximation 

A commonly used approximation to solve the flow equation for at zero external 

momenta is the local potential approximation (LPA). In this approximation, one assumes 



that the effective action has the form 

"f3 



40 



62 ] 



rf A [0] =J drj d d ~ l x jify^i + K(p)} 



(3.4) 



where p = 4>i4>i/2 and V K (p) is the effective potential. The propagator of the scalar field can 
be decomposed into its transverse {Gj) and longitudinal (Gl) components: 

G^K- q) = G t (k; q) - ^ + G l {k- q)^0 . (3.5) 

The equation for the potential, derived from Eq. (13.31) by assuming <fi to be constant, then 
reads 

W K (p) = \ J^J (S^T [ 9 « R M] W - 1) G T (n; q) + G L (n; q)} , (3.6) 

with 

Gl(K; q) = f + v^ + lvM + RM • (3 ' 8) 
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with V^(p) = dV K /dp and V"(p) = d 2 V K /dp 2 . The LPA may be viewed as the leading order 
in a systematic expansion of the effective action in powers of the derivatives of the field. 
Such an expansion has been shown to exhibit quick apparent convergence if the regulator is 



appropriately chosen [39j, |57|, |60|, l63 |. 



The choice (13.21) of the regulator leads to the following flow equation for the potential 
(see App. [£} 

o K V K (p) = Kd-iK ^ (TV — 1) ^ + j, (3.9) 

where we have defined 



I = v / K / (p) + 2pK / '(p) + ^ 2 - (3-10) 



Here, n(go) = (exp(go/^) — 1) 1 is the bosonic distribution function and Kd = Sd/ (d(2n) d ) 
with S d = 2-K d / 2 /T(d/2). In the cases that we shall study below, either the second (N = 1 
scalar theory) or the first term (large iV limit) in the braces will contribute to the flow. Note 
that, because of the one-loop structure of the flow equation, the r.h.s. of Eq. (13.91) naturally 
separates into a "thermal contribution", which involves the statistical factors and which 
vanishes when T — > 0, and a "vacuum contribution" that contains no statistical factors. 



B. Truncation and elementary analysis of the flow equations 

In the next section, we shall present numerical solutions of the flow equation, Eq. (13.91) . 
However, much insight can be gained by considering a simplified version of this equation, 
obtained by expanding the potential V(p) around p = so that m 2 = V^(p)| p=0 and g 2 = 
V"(p)\p=o/S. We shall give only the formulae for the case N — 1, but they easily generalize 
to arbitrary N. We start from the flow equation (13. 9ft and insert a truncated potential of 
the form 

V K (cf>) = K + ^<p 2 + gl<p* + h 2 js + ... 

= V K + mlp + 4glp 2 + 8h 2 K p 3 + ■■■ (3-11) 

Truncating the series at order 0(p 3 ), and neglecting h K and higher order coefficients, result 
in replacing V^(p) + 2pV tv ''(p) + k 2 by k 2 + + 2Ag 2 K p. We obtain the following set of coupled 
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Figure 5: Diagrammatic illustration of the r.h.s. of the flow equation of the effective action, 
Eq. (|3.12p . The crossed circle represents an insertion of the regulator d K R K) and the thick line a 
full propagator with mass m K . 




Figure 6: Diagrammatic illustration of the r.h.s. of the flow equation for the 2-point function, 
e.g. the mass, Eq. f)3. 13|) . The black dot denotes the four-point function (the coupling g K ) and the 
thick line the full propagator G K (including the regulator, see Eqs. (|3.7p . (|3.8p ). The circled cross 
represents the insertion of d K R K . 

differential equations at orders p°, p, and p 2 : 

0(p°): d K V K = ^_^ 1 + 2n(eK) , (3.12) 
„ ; h „ o .1 + In (e K ) — 2e K n'{e K ) 

O p 1 : d K m 2 K = -6g 2 K K d ^K d ^ , 3.13 

K 

0{p : d K g K = 27g K K d „ 1 K a 3.14 

with the short-hand notation e K := \J k 2 + m 2 K . The right hand side of Eq. (13.141) would 
change if one wanted to include the coupling h K (of order g®), but Eqs. (13. 12[) and (I3.13P 
would not. These equations have a simple interpretation in terms of the diagrams depicted 
in Figs. El and [7J these diagrams can be easily calculated with the formulae given in 
Appendix A for the finite temperature loop integrals. 

The general character of the flow can easily be inferred from a simple analysis of 
Eqs. (I3.12p . (13.131) and (13.141) . For instance, the coupling follows a four- dimensional flow 
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Figure 7: Diagrammatic illustration of the r.h.s. of the flow equation for the 4-point function, 
showing the contribution to the flow of the coupling g K included in Eq. (|3.14p (the contribution of 
the 6-point function is ignored in this equation). Black disks represent the coupling g K . The three 
channels are included in (|3. 14[) by the proper factor 3 (external momenta vanish). The crossed 
circle represents an insertion of d K R K , and the thick line a full propagator. 

from A to k ~ 2ttT, then a three dimensional flow before it freezes when k ~ m^. This 
can be seen from Eq. (13.141) : At large values of k, the bosonic distribution function can be 
neglected and e K « k so that d K g 2 oc g\j 'k and we recover the characteristic logarithmic 
flow in d = 4 dimensions. In the opposite limit k — > 0, the value of e K is limited by the 
thermal mass e K ~ mo ~ gT (see Eq. (12 .7p ). and the flow becomes d K g 2 K oc g^n 4 , leading 
quickly to a constant value for g 2 . In the intermediate range when we still have k > m K and 
thus e K m k, but one can already expand the thermal distribution function for k < T (see 
Eq. (I5.12p below), one finds d K g 2 K oc g^/n 2 , which is compatible with a three dimensional 
flow, with g K ~ k at small k. This qualitative behavior will be confirmed by the numerical 
results presented in the next section. 



IV. NUMERICAL RESULTS 



The numerical integration of the flow equation (13. 9p (using the Runge-Kutta method 
with adaptive step-size) is hampered by the flow of the potential over several orders of 
magnitude. In order to keep the zero-temperature flow under control, one can introduce 



dimension-less variables 



401 ]. but these will not be suitable for the thermal flow, which 



freezes in dimensionful variables (but would diverge in dimensionless variables as k — ► 0). 
However, one can take advantage of the fact that a large part of the flow cancels between 
thermal and vacuum contributions, as will be described in subsection IIV C\ and obtain useful 
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Figure 8: Comparison of the flow of the coupling at zero temperature for g = g(2irT) = 0.0994 
(9 A = 0.1) and for g = 1.1032 (g\ = 4). At weak coupling one observes that the 2PI /3-function 
differs form the one-loop result by a factor 1/3. At larger coupling the Landau pole is observed in 
the running of the coupling. 
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Figure 9: Same as Fig. [8] for g = 0.699 (c?a = !)■ 



results for the potential, and therefore the pressure, mass, and coupling. In the first part of 
this section, we present results for the scalar theory with N — 1 component, and in the last 



subsection for the large N limit where the LPA becomes exact 64J, l65 |. 



A. Coupling 



Figures M and [9] show the LPA flow of the coupling g K at zero temperature, for different 
values of the g = g(2irT), comparing with one-loop results as well as with those of the self- 
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consistent approximation described in Sec. [TTJ The three plots correspond to T = A/100, and 
the numerical correspondence between g and bare coupling g& is indicated in the caption. 
By construction, the curves cross each other at the scale k = 2nT. As shown by the left 
panel of Fig. [8] the LPA /^-function agrees with the one- loop /^-function, Eq. ( 12. 5p . at weak 
coupling, reflecting the scheme independence of its first coefficients. Also visible on Fig. [8] is 
the fact, already mentioned, that the 2PI resummation scheme leads to a /^-function that is 
smaller by a factor 3. As discussed in section Hi Bl at larger couplings (right panel of Fig. [Sj), 
all three curves exhibits the presence of the Landau pole where g K diverges. Note that the 
perturbative estimate of the position of the Landau pole in Eq. ( 12. lip is exceeded by the 
NPRG flow, and the moderate 2PI flow predicts a Landau pole that lies further away. The 
one-loop calculation in Fig. [8] (right panel) does not make sense for the large value of the 
coupling since then, as seen in the right panel of Fig. [SX Al < A. Figure [9] shows the flow at 
some intermediate coupling g(2nT) = 0.699 (#a = 1)- 

Thus, the presence of the Landau pole limits the applicability of the NPRG to not too 
large couplings, as already discussed in Sec. Ill Bt There is however another limitation that 
comes from the fact that the denominator in ( 13.91) . given by ( 13.101) . becomes imaginary. 
This is to be expected to happen when the coupling gets too large, as can be seen from the 
perturbative estimate: V^(0) = m\ rs — g\k 2 /{2-n 2 ) (see Eq. f !5.8j) below). 

In Fig. [10] we turn on the temperature, and compare the vacuum flow with the flow 
obtained at a temperature T = A/100. At weak coupling (left panel of Fig. [TOl . the scales 
gT and T are well separated, and we observe the general features discussed in Sec. IIII Bl 
From k = A down to 2irT, the vacuum (T = 0) flow of the coupling and the thermal 
flow agree. While the vacuum flow continues to decrease logarithmically all the way down 
to k = 0, below k = 2nT the thermal flow starts to deviate from the vacuum flow in 
a dimensionally reduced manner. The three-dimensional flow stops at the scale k ~ gT, 
where the thermal mass is generated. The situation does not change qualitatively for larger 
values of the coupling in the right panel of Fig. [H]] (enlarged in Fig. [TTJ) : The main difference 
is that the range of the three-dimensional flow shrinks as gT approaches T. We shall return 
to this at the end of the next subsection. 
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Figure 10: Flow of the coupling at finite temperature for g = 0.0994 (g\ = 0.1), left, and for 
g = 1.1032 (<7a = 4), right. For the smaller coupling, the value of g K drops only moderately 
from k = A to 2irT and the thermal coupling reaches a fixed value at g K =o = 0.093. For the 
larger coupling, the beta function forces a steep drop from g\ = 4 to g = 1.1032, and further to 
g K =o = 0.68 for the thermal coupling. 
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Figure 11: Enlarged view of Fig. [10] (right) to show the region of dimensionally reduced flow. 

B. Mass 

The freezing of the flow at soft scales is caused by the thermally generated mass 
m 2 K = V^(p = 0), shown in Fig. [121 As is the case for the coupling, the flow in the range 
2ttT < k < A agrees between vacuum and thermal flow. The system starts at a large neg- 
ative mass squared, which can be estimated perturbatively as m\ = —g\A 2 /(2n 2 ) + 0(g\), 
see Eq. (15 .8p . In the course of the flow, the mass acquires positive contributions through 
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Figure 12: Flow of the mass for g = 0.0994 (g A = 0.1), left, and for g = 1.1032 (g A = 4), right. 
The mass reaches the value m = 0.094T for g = 0.0994 and m = 0.622T for g = 1.1032 as k — > 0. 
(For T = 0, the quantity plotted is m 2 ,^ K /T 2 .) 
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Figure 13: Thermal mass for scalar theory with N = 1 component as a function of the coupling. 
The result obtained from the non-perturbative renormalization group labeled "RG-LPA" gets close 
to the result obtained by 2PI resummation. Both show moderate behavior where the perturbative 
contributions represented by curves of order g 2 and g 3 fail to converge. 

quantum fluctuations. The initial condition m\ has been tuned at zero temperature, such 
that vanishes. At finite temperature, the flow of the mass deviates below the scale 

k = 2nT and builds up the thermal mass m ~ gT that can be observed in Fig. [121 

The values obtained from the flow equation in Figs. [10] and [12] are combined to produce 
the result displayed in Fig. [131 For example, starting the numerical calculation from g\ = 4, 
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one obtains a coupling g = g(27iT) = 1.1032. At this value of the coupling, the mass takes 
the value m/T = 0.622 (Fig. [T2l in the physical limit k — > 0. This information is put 
into Fig. [131 an d the procedure is repeated for other values of the coupling. The curve 
thus obtained, labeled "RG-LPA", is smooth, in contrast to those displaying the results of 
strict perturbation theory through orders g 2 and g 3 . For comparison, in Fig . [TBI we also plot 
the results obtained through the 2PI resummation described in Sec. II [17] . The difference 
obtained between the RG-LPA and 2PI approaches has two origins: on the one hand, the two 
calculations correspond to two different /3-functions: /3 2 pi = Ame-ioop/3 (see Fig. [SJ). On the 
other hand, the two calculations use different renormalization schemes: the 2PI curve shows 
the scale choice in the MS-scheme with /2j^g = 2nT . In comparison, the NPRG approach 
defines the coupling from the vacuum flow at k — 2nT and the quantities at other scales 
depend on the choice of the regulator. 

From Figs. [10] to [13] we can already extract one of the main results of our work. Both 
in the weak and the strong coupling regimes, the scale dependence of the running coupling 
and of the mass exhibits a similar behavior: in both cases, thermal effects start to play a 
role around k ~ 2tcT and saturate at k ~ mjj ~ gT. When going from weak to strong 
coupling, there is a competition between two effects: on the one side, the range of the values 
of k where thermal effects are important shrinks because the thermal mass increases with 
the coupling; on the other side, the amplitude of the three-dimensional flow increases with 
increasing values of g. Both effects are visible in Figs. [lU]and[TTJ The net effect is that the 
variation of the coupling along the three dimensional flow is limited; for instance, in Fig. [TT] 
for g = 1.1032, it amounts to a 30% effect. 

One can also observe in Fig. [TTIthat. at strong coupling, the flow saturates at a value of 
k < gT. This is because the thermal mass reaches a value lower than gT, as illustrated in 
Fig. [131 In fact, one may argue that the thermal mass cannot exceed a value ~ 2irT, since 
otherwise all thermal effects would disappear, which is in contradiction with the existence 
of such thermal mass. 



C. Pressure 



Figure [H] shows the flow of the pressure P = —V(p = 0) for weak and strong couplings, 
normalized to the free pressure Pq = (it 2 /90)T 4 . We choose Va(p = 0) such that the zero- 
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Figure 14: Comparison of the flow of the pressure in the weak-coupling regime for g = 0.0994 
(d A = 0-l)> left panel, and in the strong-coupling regime for g = 1.1032 (g\ = 4), right panel, with 
T = A/100. The additive constant is chosen such that the vacuum pressure vanishes as k — * 0. The 
thermal pressure reaches the value of P = 0.9984P for g = 0.0994 and P = 0.945P for g = 1.1032 
as k — > 0. 

temperature pressure P k= q(T = 0) vanishes. This requires however a delicate cancellation. 
Indeed, a perturbative calculation (see Eq. (15. 7p below) shows that the vacuum pressure at 
k = A is of the order A 4 . Thus, for example, both curves for the g = 1.1032 plot start at 
Pa = —415, 470.518Po and acquire a relative difference during the flow from k = A to k = 
of 0.945Po- Clearly, this cancellation over six orders of magnitude poses a challenge to the 
numerics. One should note though that the thermal flow above k > 27rT lies exponentially 
close to the zero temperature flow due to the bosonic distribution function in Eq. (13.91) . 
and the bulk of the thermal contribution builds up between k ~ 27rT and k ~ gT. This 
difference between thermal and vacuum flows is shown as a light curve in the plots. It is 
therefore sufficient to start the integration of the thermal pressure at some intermediate 
scale 27rT < C A. The zero-temperature and finite-temperature flows of the potential 
only reach a maximum value of the order of A 4 A 4 and their difference is numerically 
much easier to handle. On the other hand, Ai can not be chosen too small, or contributions 
to the result will be neglected. Practically, a value of Ai = 20T turned out to be a good 
compromise between loss of accuracy due to neglecting k > Ai contributions and gain of 
accuracy due to reduced numerical cancellation. This procedure also improves the numerical 
accuracy of the thermal mass, shown in Fig. [121 although the situation is not as critical there 
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Figure 15: Pressure as a function of the coupling (left panel) and as a function of the mass (right 
panel). As two physical, scale-independent quantities are compared, the right plot is renormal- 
ization scale independent and shows an improved agreement between RG-LPA and 2PI results. 



(the mass only grows as k 2 , as opposed to the pressure that grows as k 4 ). 

Since the pressure, mass, and coupling constant all approach constant values as k —>■ 
0, dimensionless variables that one may introduce for the zero-temperature flow cease to 
be useful: For example, as the dimensionful potential V^_>o approaches a constant, the 
dimensionless variable diverges as f K ^o ~ 1/k . Numerically, we have still followed the 
approach of dimensionless variables also for thermal quantities, as k does not need to become 
too small to measure the limiting value reached when k — > 0. 

The left panel of Fig. [15] shows the pressure as a function of g. It follows combining 
results from Figs. [TD] and [HI in the same way as Figs. [TU] and [£2] lead to Fig. [TSJ Due to 
the limitations of the maximum value of the coupling that we can use (see discussion on 
page H7|) . it is necessary to reduce the ratio T/A for the larger couplings in Fig. [15j The 
range depicted (up to g ~ 1.5) is about the range in which both pressure and thermal mass 
are still independent of this choice of T/A, up to plot accuracy. For the largest g ~ 1.5 we 
have chosen T/A = 1/20; this agrees for smaller g with T/A = 1/50 or 1/100. The value 
T/A = 1/10 could extend the NPRG method to even larger g, but then the result visibly 
starts to deviate from the T/A = 1/20 result. Similar to the plot for the mass in Fig. [THl 
we observe an improved result "RG-LPA" over perturbation theory through orders g 2 and 
g 3 . We obtain a deviation of similar magnitude when compared to the 2PI resummation 
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scheme, part of which is caused as mentioned above by different renormalization schemes. 
Note that the 2PI curve lies below the RG-LPA curve. 

We can circumvent the problem of renormalization scale by comparing only physical 
quantities. Figure [15] (right panel) displays the pressure as a function of the mass. For the 
LPA curve, both quantities, pressure and mass, are extracted at k = 0, so there is no more 
reference in this plot to a particular choice of a scale for the coupling g. The perturbative 
contributions g 2 and g 3 on the other hand can be considered as scale-independent only up 
to the order they have been calculated to. These curves may change when extracted at a 
different scale. In this plot we observe improved agreement between the RG-LPA curve and 
the 2PI curve, but also note that they do not completely agree. Particularly, at larger m/T, 
the LPA curve starts to bend up again. This is the region where the NPRG equations get 
sensitive to the Landau pole and may also fail because of imaginary parts in the denominator 
of the flow equation as discussed on page [T71 so this region should be considered with care 
(the fact that the 2PI curve does not show the same bending as the RG-LPA one presumably 
reflects the fact that the 2PI Landau pole is artificially pushed to higher values because of the 
factor 1/3 in the 2PI /5-function). The break-down of perturbation theory can be pinpointed 
where the g 3 contribution deviates visibly. This curve bends backwards, because m(g) is 
not a monotonous function at this order. Surprisingly, the g 2 curve seems to remain a good 
approximation to the non-perturbative approaches. 



D. Large N limit 

In the large N limit, only the transverse propagator in Eq. (I3.9P contributes to the flow. 
The various quantities scale in the following way as iV — > oo: 

P,V,p ~ N, 
V^p),m 2 ~ 1, 

g 2 ~ 1/N. (4.1) 

It is then advantageous to plot results as a function of an effective coupling g c g = gy/N. 
Note that both P/Pq and m/T are iV-independent quantities. 

In Fig. [16] we plot the predictions of the RG-LPA for the mass and the pressure as a 
function of the coupling constant, together with the exact large iV-limit results, which are 
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Figure 16: Comparison of mass (left panel) and pressure (right panel) at large N, normalized 
to temperature and the free pressure, respectively. The coupling g denotes the effective coupling 
9eS = gV^N in the limit N — » oo. 
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Figure 17: Pressure as a function of the mass in the large ./V limit. In the comparison of scale- 
independent quantities, the RG-LPA numerical result exactly reproduces the large N limit. 

known (see e.g. [l|). In the large N limit, the LPA approximation also gets exact 64, 65]. 
The reason why the curves of Fig. [16] do not exactly agree is because of the already mentioned 
scheme dependence. If one plots only scale-independent quantities, like the pressure as a 
function of the thermal mass, then one recovers indeed a full agreement between the RG-LPA 
and the known exact large N limit results (see Fig. [T7l) . 
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V. PERTURB ATIVE EXPANSION OF THE FLOW EQUATION 

In the spirit of other resummation schemes proposed in the literature, we show that 
the NPRG in the LPA, reproduces perturbation theory up to order g 3 . Of course, the 
full numerical evaluation of the flow equation in Sec. IIVI includes higher orders in a non- 
perturbative way. 

We start from the truncated flow equations (I3.12p to (13.141) where the neglected term in 
(13. lip . h K , is of order g & K . These equations for the mass (13.131) and coupling constant (13. 14j) 
from a closed set of equations that can be solved order by order in g. The pressure can then 
be calculated using the results for m K and g K . 

In the case of the vacuum flow, the perturbative hierarchy is quite clear: m 2 , is always 
~ g 2 K 2 , and one can therefore expand e K = \J k 2 + m 2 ~ k for all k. 

At finite temperature, on the other hand, a thermal mass of the order gT is generated, 
and expanding for small m K necessarily fails when n < gT. In this case, it is useful to 
introduce an intermediate scale £ that lies between soft and hard momenta gT < ( < T 
(e.g. £ ~ ^JgT). Different expansions can then be performed in these two regions: For hard 
momenta, n > £, 

el = K 2 + m 2 K = K 2 (l+O((gT/0 2 )). (5.1) 
For soft momenta, k < £, we can expand the thermal distribution function 

nU) = --- + — + ZoU^/Tf). (5.2) 
K J k 2 12T k VVW ; ; V ; 

The leading term, T/k, trades one dimension of a (i-dimensional integration by the temper- 
ature, thereby leading to a dimensionally reduced system. From the calculation of section 
III A| one expects order g 2 contributions to arise from integrals dominated by hard momenta, 
and g 3 effects to build up at soft momenta. However, as we shall see, delicate cancellations 
occur in the course of the calculation. In order to cleanly reproduce strict perturbative 
results at order g 3 , we shall expand the actual running quantities g K or m K around their 
leading order perturbative contributions: 

2 2 , r 2 

9k = 9 +og K , 

m 2 . = m 2 + 5m 2 K , (5.3) 

with m 2 = g 2 T 2 (see ( I5.22p ). where g without subscript is to be understood, as before, as 
9 = 5 , vac('« = 27rT). As we shall see, in the region k > £ the correction 5g 2 K will be formally of 
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order g 4 (see (15. lip ), while for k < £ the correction Sg 2 will be formally of order g 3 (see ( IB II) ) 
and 5m 2 . of the order g 2 ^T (see ( IB 2 ft and f)B3[) ) . Technical details that amount to 0(g A ) or 
0(g A In g) contributions are deferred to App. [Bj the LPA is not expected to reproduce them 
correctly. 



A. Zero-temperature flow 

At zero temperature, all the distribution functions n(e K ) and their derivatives vanish in 
Eqs. ( 13.121) to ( I3.14p . and we are left with the following set of differential equations, that we 
expand in powers of g 2 , 

^^^-^(^-^J+OWS, (5.4) 
d.ml = - 6g2 " K *- lKd = -6glK d ^ + 0(g*) , (5.5) 



27g i K K d . 1 K d 

c5 



« = -^-P— = 21glK d ^K^ + 0(g») . (5.6) 



In these equations, we anticipated that the mass m 2 will be of order g 2 (i.e. we assume that 
m K=0 vanishes, so that m 2 is entirely due to quantum fluctuations). For d — 4 the solutions 
of these equations read: 

V» = V A + ^4- + 0(g 2 A ), (5.7) 



48vr 2 



K 



2 



< = -fcf + ofol), ( 5 - 8 ) 



: 2tt 2 

In deriving these equations, we have assumed that jaCI, and also taken into account the 
possible occurrence of large logarithms in getting Eq. (I5.9p . To get Eq. (15.81) . one integrates 
the r.h.s. of Eq. (15.51) by parts using (15.61) . Finally, the integration constant has been adjusted 
so that m 2 =0 vanishes, i.e. we have m\ = —g\A 2 /(2ii 2 ). Using these results, we can actually 
proceed one order further in the potential: 

V K = V A + + 9 ^ ~ ^ + 0(g{). (5.10) 

48vr 2 192vr 4 vyA; v ; 
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B. Coupling at finite temperature 

In order to find the perturbative solution of Eq. (I3.14p we introduce an intermediate scale 
£, as discussed above, with gT < ( C T. For k > £ thermal effects are subleading and g 2 . 
can be obtained directly from (15. 9ft 

^ = » 2 + ^l°g^ + 0( 9 «). (5.11) 

In the region k < £ we can expand the distribution function in Eq. (13.141) according to (15. 2p 
which gives a remarkably stable expansion (i.e. all powers from e° K to e K drop out) 

2 1 fiT f 6 

1 + 2n (e K ) - 2e K n'(e K ) + -e 2 K n" (e K ) = — + + 0(p . (5.12) 

The flow equation then becomes 

« = (5-13) 

with e 2 K = k 2 + m 2 . By writing m 2 = m 2 + 5m 2 (see Eq. (15. 3p ). and expanding out the 5m 2 K 
correction, we obtain for the leading contribution 

5 



9l 



3T / 3 k 5k + 3/«m z 
■ arctan 



7r 2 \m m (k 2 + m 2 ) 2 



+ 0(£/(gT)). (5.14) 



In the limit k' — > and gT <C ^ — > oo and expanding the result in g, we recover the result 
obtained earlier 

dLo = 9 2 - f <7 3 + Otft/T) . (5.15) 

(The expression for general k is given in (IBlh ) . Note that the result is, to order g 3 , inde- 
pendent of the temperature T. This may seem at first sight surprising. However, the only 
possible effect of of the temperature is a change of the renormalization scale, and this only 
influences the result at order 0(g 4 ). A careful calculation through order g 4 (presented in 
Appendix IB lj) reveals how the value of the coupling at k = is connected to the coupling 



at k > T: 

2 2 9c/3 9g* ( k 61 32 . / 4 - 

£=° = £-^-^{lo e — +J -- + - 2 ) + 0( A ), (5.16) 

where 7 is Euler's constant. Note that this result is obtained from perturbatively expanding 
the flow equation within the LPA. 

Figure [18] shows a comparison with the numerical solutions. The curves labeled "g 3 " 
show the sum of the g 2 value at k — 2ttT and the K-dependent correction of order g 3 from 
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Figure 18: Same as Fig. [TUl compared to analytic solutions. Dashed lines give perturbative solution 
through order g 3 (flowing with k), with g taken at k = 2ttT, and g 4 (only go value), with g taken 
at k = 2ttT. For the left panel, also the solution of truncated LPA is given by the light curve. 

Eq. (1B1I) . This formula is strictly speaking only valid in a region k < £ and introduces an 
error of the order 0(gT/£). This error is the reason why the g 3 curve does not agree with 
the numerical value of g at k = 2ttT, but reaches this value only for k — > oo. Also, through 
order g 3 there is no renormalization scale dependence, and the analytical solution does not 
follow the vacuum flow of the numerical result. The curves labeled "g 4 " in Fig. [TBI show the 
result (I5.16P with g(n = 2irT). 

For the left plot of Fig. [18] also the result of truncated LPA, obtained by numerically 
integrating Eqs. (13. 12[) to (13.141) . is given. While for small coupling in the left plot of Fig. [18] 
the truncated LPA misses the full numerical result of the scattering amplitude g K =o by 
0.3%, at the larger coupling g(A) = 4, the truncated LPA would underestimate the value 
of g(n = 2ttT) by 28% (not shown in the plot). If we instead lower the cutoff A such that 
the vacuum couplings at k = 2nT agree, the truncated LPA will give a result for g K=Q that 
is 23% larger than the full numerical result. Thus, while this truncation provides reliable 
estimates at small coupling, it cannot be used at strong coupling. 
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C. Thermal mass 

We consider now finite temperature effects and focus first on the thermal mass. This is 
defined by separating the mass in equation (13 . 1 3[) into vacuum and thermal pieces: 

m vac,K = "£(T = 0), (5.17) 

m th,K = m l( T )- m l aC ,K- ( 5 - 18 ) 

The thermal contribution can then be calculated using e R : = y / ' k 2 + m 2 hK + ml acK , and 
reads: 

™th,«' = -6^-i / dKK \ g K - 3 ^ac,K , K 2 7 ^2 f ' (5 " 19 ^ 

1. Thermal mass to leading order 

As we have seen earlier, the g 2 contribution to the thermal mass is dominated by hard 
momenta. Using Eqs. (15. ip and (15.31) . we obtain as leading contribution 

,A 

m t 2 hjK , = l2g 2 K d ^ / dn (n(«) - «n'(«)) + 0{g 4 ) . (5.20) 

in' 

Although, strictly speaking, this expression is valid only for k' > £, it is IR safe, and the 
error introduced by sending n' to is of the order 0(g 2 ^T). On the UV side, the error 
introduced by replacing A by oo is exponentially small (it is of the order 0(TAe _A//T ) for 
the part with n, and 0(A 2 e _A//T ) for n'). We finally obtain for k' = 

POO 

m t 2 h = 12g 2 K d _ 1 (d - 1) / dn K d - 3 n(^) + 0(g\T) + 0(AV A/T ) . (5.21) 

Jo 

For d = 4 dimensions, this gives 

"4,o = ™ 2 = S 2 ^ 2 • (5-22) 

Thermal mass at next-to-leading order 

In order to obtain the next-to- leading order in Eq. (I5.19p . we need to treat the small 
k limit carefully. The resulting contribution, of order g 3 , is calculated by subtracting the 
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leading g 2 contribution of Eq. (15.201) from the unexpanded Eq. (15.191) . 



'th,0| 9 3 




2 1 + 2n(e K ) - 2e K n'(e K ) 




(k 2 + m2 acK )3/2 



1 



(5.23) 



Note that we keep the ^-dependence in the first line, as the flow of g K is affected by g 3 cor- 
rections (see Eq. (15.1 5ft ) that may change the mass at order g 3 . In fact, a careful calculation 
reveals that such contributions would affect the result only at order g 4 In g (see IB 2j) . 

Since the dominant contribution to the integral (15.231) stems from the infrared region, 
we split the integration range by an intermediate scale £ with gT < ( < T. For hard 
contributions n > £ one can expand m K , m racA (< k) in the denominators and in e 2 K = 
K 2 + m 2 , and use g^ ac K = g 2 + 0(g 4 ); one observes then that the result vanishes up to 0(g A ). 

For soft contributions < k < £ we can expand the thermal distribution function (15.21) . 
To extract the g 3 term, it is sufficient to keep the leading term of (15.21) and use coupling and 
mass only up to order g 2 , i.e. g K — > g and m K — > m = gT. Equation (I5.23P then simplifies to 



For the thermal pressure we expand equation (13. 12ft in a similar way, subtracting the 
vacuum piece 




(5.24) 



We recover the correct result for the correction to the mass at order g 3 . 
We can finally write the perturbative thermal mass as 



ml,o = T 2 9 2 ~-9 3 + 0(g 4 logg) 



(5.25) 



D. Pressure at finite temperature 



Pth = -Vth.0 = -(K(T) - V K (T = 0))U =0 , 



(5.26) 



and obtain from Eq. (13. 1 2[) 
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where as before e 2 K := k 2 + m 2 h K + m^ ac K = k 2 + m^. The leading contributions at order g° 
and g 2 come from the region k > £, in which both, m t 2 h K and m 2 ac K , can be expanded from 
e K (see (EH])). We get 

,A 

Pth = y ^{^"^(k) (5.28) 



_d-3 



Kk + 2 (m t 2 hjK + mLJ (n(«) - kn'(K))] (5.29) 
+ O((^T/0 4 )} (5.30) 

where the first line is of order g and the second line is of order g 2 . The third line in this 
expansion scheme is formally of order g 4 , but we anticipate that soft modes with k < £ will 
contribute a term of order g 3 . The expressions (15.281) and (15.291) are IR safe, and £ can be 
sent to in these expressions without affecting the perturbative result through order g 2 . 



1. Free pressure 

The free pressure is easily evaluated from the first line of Eq. (15.281) 



p = lim Kd-i / dKK d 1 n(K) 



T -^U d {l)T d (5.31) 



7T 

where Li^(z) is the polylogarithm function. This gives the known values for the free pressure 
in 2, 3, or 4 dimensions 

P-,3. = ^,^,^. (5.32) 

2. Pressure at next-to-leading order 

The contribution to order g 2 is already harder to obtain. This is because, as Eq. (I5.29j) 
shows, both, m 2 h K and K , contribute on equal basis to the g 2 result of the pressure. Of 
course, at the end of the flow, the zero-temperature mass vanishes m 2 AC K = 0, as we have 
tuned the mass m\ just in this way, but this does not mean that we can neglect the effect 
of m 2 &c K during the evolution of the potential from n = A to 0, even if the final result turns 
out to be strictly proportional to the thermal mass p\ g 2 oc m 2 h . 
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We shall evaluate Eq. f!5.29j) in d — 4 dimensions in three pieces Pth|„2 = P^*™^ + 



(vacxth) . (thxth) .,, 
Pth +Pth Wlth 



(thxvac) [ j 2 (r- oo\ 

Pth = — T / dKKm th, K , (5-33) 







4 

^(vacxth) = _3 j d/cKTT^^/c)-^^)), (5.34) 

^(thxth) = -^1 J dKKm 2 thK (n(K)-kn'(K)). (5.35) 
For the first term one uses the explicit leading-order expression (I5.20p for m 2 h : 

t>A pA 

Pth hXvac ^ = — 3g 2 K 2 / dtz k j dn' k' (n(K') — k' n' {k')) . (5.36) 



o 



One can integrate this expression by parts, obtaining for A — > oo 

p (thxvac) = _3g^Ki f™ dK (^ n ( K )_ K ± n >( K )} = _^K! f°° dKK z n ( K ) 

g 2 T* 
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This piece is exactly canceled by the following piece of p t h 



(5.37) 



3 (vacxth) = J_ I dKKm 2^ n ^ _ Kn '( R )) = +— — , (5.38) 



vacv V / V J / 1 ^2 ' 

where we have used Eq. (15. 8p for the mass m 2 ac = —g 2 K 2 /(2ir 2 ), again in the limit A — > oo. 

This cancellation reflects a general property. In a perturbative calculation, the terms 
mixing vacuum and thermal parts may contain ultraviolet divergences coming from vacuum 
subdiagrams. Since there can be no di verg ent terms involving the temperature in the final 



results, such terms must cancel 
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66, 
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The residual contribution to p t h therefore solely comes from 

pfjfcxih) = j dKKm 2 th (n(K) - Kn'(K)) . (5.39) 
Identifying parts of the integrand with the flow equation (13. 13[) for the thermal mass 

d K m 2 th = -l2g 2 K 3 K(n(K) - Kn'(n)) , (5.40) 

we can write the final contribution as 

th I f A 1 g 2 T A 

Pth| g2 = p^ xth) = ^2 J dKm thd,m 2 h = -^"4 = — ^g- • ( 5 - 41 ) 
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3. Plasmon contribution to the pressure 

In order to calculate the g 3 contribution, we proceed analogously to the calculation of the 
mass at order g 3 : We form an IR dominated integral by subtracting the known contributions 
at order g° and g 2 in Eqs. (15. 28ft and (15. 29ft from the exact integral Eq. ( 15.27ft : 

1 + 2n(e K ) 1 



Pth\ 



93 \L " \" I 2e K 2J^+mi 







K d -i J dn < 


H 






-K d ~ l n{K) + 


K d-3 


4 



This integral is dominated by small values < k < £ and we can replace the upper integra- 
tion limit A by £ without changing the result at order g 3 . We proceed as before, expanding 
the distribution function as in Eq. (j5.2p and neglecting m^ acK ~ 0(g 2 n 2 ), and obtain 



2 , 2 



Pth\ g3 = TK 3 I oIk^ — -k + m thK 



TK^ I dn— = TKj, m 3 arctan — 

K + m 2 m 



(5.43) 



where in the second line we have expanded the mass as in (15.31) . Sending the upper inte- 
gration limit £ to 00 in (I5.43P changes the result only beyond order g 3 , and we can finally 
write 

m 3 T T 4 

PthLa = = 9 3 • (5.44) 



VI. CONCLUSIONS 



In this paper we have applied the non-perturbative renormalization group and its local 
potential approximation to a scalar field theory at finite temperature to calculate its pressure, 
screening mass, and scattering amplitude, covering both weak and strong couplings. In the 
perturbative regime, we have shown that the LPA reproduces perturbative expansions of 
pressure and screening mass up to and including the order g 3 . The latter is important as 
thermal effects at weak coupling are dominated by the effects of order g 3 , and we have 
indeed verified that the three-dimensional flow of the coupling constant is entirely given by 
these effects at small couplings. The LPA allows a smooth extrapolation of weak coupling 
results into the regime of strong coupling. Of course, because of the presence of the Landau 
pole, arbitrarily large values of the coupling cannot be reached. However the range of values 
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of g that can be explored allowed us to demonstrate a clear improvement over the strict 
expansion in terms of the coupling constant. 

We have compared the results obtained within the LPA with those of a simple 2PI 
resummation: both methods lead to very similar results in the extrapolation to strong 
coupling. This is not too surprising if one considers the diagrammatic content of the two 
approximations, and the fact that both become exact in the large N limit. The LPA has 
the advantage over the 2PI formalism that it yields the one-loop beta function correctly 
(while it is necessary to go to higher order to achieve that in the 2PI formalism). This is 
because the three channels of the 4-point function are treated simultaneously in the LPA, 
albeit approximately (loop insertions carry no external momentum in the LPA). 

The NPRG provides an understanding of what happens as we move to strong coupling. 
As we have seen two effects compete as g increases: the region of three dimensional flow 
shrinks because the thermal mass increases, and the amplitude of the three dimensional 
flow grows. In a sense dimensional reduction continues to play an important role, but is no 
longer related to a weak coupling effective theory. It is the competition between these two 
effects that is responsible, within the LPA, for the stability of the results obtained for the 
pressure with increasing couplings: the corrections to the pressure remain modest; even for 
the largest available couplings they never exceed 10%. 

Of course, the calculations presented in this paper have limitations. We have already 
mentioned the impossibility to go to too large couplings in scalar field theory. Such a 
limitation would not appear in an asymptotically free theory like QCD where many of these 
difficulties of the perturbative expansion that we have discussed are encountered (among 
others). Another limitation comes from the fact that the LPA ignores the momentum 
dependence of the self-energy, and the effect of the width of quasiparticles is therefore 
entirely neglected. Treating such effects is beyond the reach of the derivative expansion. 
However the technology to extend the present study in this direction exists 
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Figure 19: Integration path for Eq. (|Aip , 
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ence on the Exact Renormalisation Group, in Lefkada, Greece, September 18, 2006, D. Litim 
and J. Pawlowski informed us about their own investigation on this subject |74j . 

Appendices 

Appendix A: FLOW EQUATION AT FINITE TEMPERATURE 

Sum-integrals that occur in various parts of our work can be calculated from standard 
contour integrations using the formula: 

(Al) 

with the Matsubara frequencies u n = 2nnT, n(q ) = (exp(q /T) — 1) _1 and n(q ) = 
(exp(q /T) — 1) _1 the bosonic distribution function. Equation (lAlj) is valid if f(qo,q), 
considered as a function of go has singularities only on the real axis and is an even function, 
/(<?()) q) — f(~Qo, q)- The contour C\ is displayed in Fig. [19j The following useful formulae 
are easily established: 

T y 1 = 1 + 2n (^) (A2 ) 

TV L_ = l + 2n«-2o;X 
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TV 5 = -^—Sl + 2n K ~2u K n'+-u 2 n''\. (A4) 

As an application, let us derive Eq. (13.91) for the effective potential. Consider the longitu- 
dinal modes first. With the regulator (13. 2p the q integration in Eq. (I3.6P can be performed 
analytically, leaving 

d K V K (p) = K d ^K d Tj2 5—1777 u 1 , T/ „, ri 2 > ( A5 ) 

where K d is defined after Eq. fl3TTUjl . When Vj.(p) + 2pV"(p) + k 2 > 0, we can use 0X2]) and 
obtain 

<9 K K(p) = ^-iK d 2n( ^ + 1 , (A6) 

where ui^ has been defined in ( 13.101) . The transverse piece that leads to ( 13. 91) is calculated 
analogously. 

Appendix B: ANNOTATIONS TO PERTURB ATIVE CALCULATIONS 

In this Appendix, we collect some technical details of the perturbative expansion of the 
flow equations. We will use the following results through order g 3 : 

The coupling at finite temperature for arbitrary soft k < £ follows from (15. 14j) and reads 
at order g 3 : 

£U = % farctan^ - I - gg^ + V^ + + O (( ff T/0^ . (Bl) 

The thermal mass has the following behavior for small k which can be obtained by 

expanding (15.201) and evaluating (15.241) for arbitrary k 

ATk k 2 
■' V ^ 2 ~ + 2^ 2 



<*L = 9 2 [T 2 -^ + ^ + 0^)), (B2) 



= ^ 1 " + f + °(^) 2 ) ) • CB3) 

The vacuum flow m^ ac K is given in ( 15.81) . 

1. Coupling at finite temperature through order g 4 



In order to calculate the constant under the logarithm of the g 4 coefficient for the thermal 
mass, we have to be more careful in integrating the flow equation for the coupling (13.141) . 
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Dividing this equation by g\ and multiplying by dn, we can integrate the left and right hand 
sides of ( I3.14p independently. The r.h.s. contains IR as well as UV problems that we have to 
cope with. We introduce an intermediate scale gT < ^ < T (e.g. £ ~ ^/gT) and calculate 
this integral in two pieces. For k > £, we can neglect m K oc gT <C £ in e K and regulate 
the remaining integral by subtracting its leading divergent pieces in the IR and UV limits. 
Abbreviating X(e K ) = (1 + 2n(e K ) - 2e R n'(e K ) + jefri" {e K )) / (%., we get 



9 r , 

dK 



2tt 2 



k 4 X(k) - ^-9(T - «) - - T) 

3 fir k 



2vr 



9 (4 + 7-log(27r)) (B4) 



with Euler's constant 7 = 0.577216. This result is obtained by integrating those pieces of 
X(k) which contain derivatives of n(n) by parts and using for n(n)/ k 

ro ° dK ( 2 2 + k \ . . 

which can be derived by multiplying the terms by K a and analytically continuing the result 
to a — » 0. Using ( 1B4I) . we can write the «; > £ contribution as 

JL jT* «fc = A (HI + log ^ + 7 - I + Ott/T)) + 0(A (B6) 

valid for sufficiently small ^ C T and sufficiently large A ^> T. 

For k < £, we can expand A(e K ) as in ( 15.121) . In a first step, we use the expression (15.14ft 
that leads to the g 3 contribution of the mass and expand it for large £ ^> gT 

24T ft , k a 9T 24T T 



, , dn 5 = - + — O((m/0 3 ). (B7) 

It is assuring that the l/£ terms cancel when adding (1B6|) and (1B7|) . but we are still missing 
an 0(1) contribution which enters on the one hand through the missing k dependence of 
m in ( 1B7I) . and on the other hand through the missing g 3 contribution to m. Both can be 
taken into account by replacing m 2 — > m 2 = m 2 acK + mj^ J^ 3 using (15. 8ft . ( 1B2I) . and (1B3I) . 
The missing contribution can be thus calculated by 

24T Z^ 00 , r k 4 
hm — — / dK 



g^O 7T 2 



k 2 + m 2 ) 3 (k 2 + g 2 T 2 ) 3 _ 

— / d« f + T (T + k ) ( — h — — arctan - ) + O j 

7T 4 7 (t 2 + k 2 ) 5 V 3 V V2 T K T J JJ yy> 

48 97 

7T 4 87T Z 
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where we have expanded the first denominator around k 2 + g 2 T 2 and substituted R = tt/g 
to form the g — > limit. 

We obtain the final result by combining (1B6I) + (1B7I) + (1B8H : 

-At f A dKK A X{e K ) = l-(* + log— + 7 - — + ^)+0(g) (B9) 
2vr 2 7 V ; 2tt 2 V £ 2vrT ' 12 3tt 2 / W V ; 

Inserting this result into (I3.14p . we obtain 

11 9 9 A A 7 32 \ . . 

5? - 51 = + 2^ ('° g 2^ + -< - 12 + SPj + °<*>- < M °> 

or, expanding and using (I2.10p to change from g& to g K for T < k < A, we obtain 

2 2 9 2« 9 9t A K 61 32 \ ^. 5 , . iN 

which connects the value of the coupling at k = to the value of the coupling k > T. 

2. Thermal mass through order g 4 In g 

When deriving the thermal mass at order g 3 , in (I5.24p we had neglected the flow of g K 
and m K at order g 3 which may change the result of the mass at order g 3 . We will see that 
in the IR region k < £ (with gT < ^ < T) such corrections would only affect the result at 
order g A lng, thereby justifying the IR calculation (15.241) . In the UV region k > £ on the 
other hand, these corrections would affect the result only at order g 4 , thereby establishing 
(I5.24p not only as the correct result for the IR, but for all k at order g 3 . 

The g 3 correction to the coupling only affects g K in the first line of Eq. (I5.23p . We can 
therefore expand g\ — > g 2 + 5g 2 with Sg 2 = g 2 \ g 3 from ( IBip . (The g 2 in the second line of 
(I5.23P is not affected by this, as we subtract a clean g 2 contribution.) In the region k < £, 
we can expand the thermal distribution function in (15.231) . The correction to (I5.24p is then 
given by 

4T , r 2 96# 4 T 2 / gT 25 9tt 2 37tgT m(M2 A .„ , 

This contribution is of the order g 4 lng. We expect the log(gT/£) from the IR region k < £ 
to combine with a corresponding log(£/T) in the UV region k > £, but this calculation 
is beyond the scope of this paper. Since g 3 corrections to the coupling and mass are IR 
results stemming from k < £, it is clear that no further g 3 contribution can build up in the 
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UV. Therefore, there can not be any additional correction to the thermal mass, and (I5.24p 
already captures the full g 3 term. 
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